Intensive thermodynamic parameters in nonequilibrium systems 
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Considering a broad class of steady-state nonequilibrium systems for which some additive quan- 
tities are conserved by the dynamics, we introduce from a statistical approach intensive thermody- 
namic parameters (ITPs) conjugated to the conserved quantities. This definition does not require 
any detailed balance relation to be fulfilled. Rather, the system has to satisfy a general additivity 
property, which holds in most of the models usually considered in the literature, including those 
described by a matrix product ansatz with finite matrices. The main property of these ITPs is to 
take equal values in two subsystems, making them a powerful tool to describe nonequilibrium phase 
coexistence, as illustrated on different models. We finally discuss the issue of the equalization of 
ITPs when two different systems are put into contact. This issue is closely related to the possibility 
of measuring the ITPs using a small auxiliary system, in the same way as temperature is measured 
with a thermometer, and points at one of the major difficulties of nonequilibrium statistical me- 
chanics. In addition, an efficient alternative determination, based on the measure of fluctuations, is 
also proposed and illustrated. 

PACS numbers: 05.20.-y, 05.70.Ln, 05.10.Cc 



I. INTRODUCTION 



The general concept of intensive thermodynamic pa- 
rameters plays a crucial role in equilibrium statistical 
mechanics. For systems in contact, ITPs like tempera- 
ture, pressure or chemical potential equalize their values 
once the equilibrium state is reached, provided that their 
associated quantities, energy, volume or number of par- 
ticles, can be exchanged. Since this property holds even 
in the case of systems that exhibit different microscopic 
dynamics, this equalization became the key criterion in 
equilibrium statistical mechanics to study the influence of 
the environment on a given system, for example when a 
reservoir is connected to it. Moreover the theory of phase 
coexistence as well as the measurement of, for instance, 
temperature with a thermometer draw on this powerful 
concept. 

Indeed this potency of the ITP formalism motivates 
the endeavor to generalize the notion of ITPs to nonequi- 
librium systems. There exist several different approaches 
mostly focusing on the generalization of temperature 
out of equilibrium. For stationary nonequilibrium sys- 
tems that fulfill a local equilibrium condition, equilibrium 
properties are locally recovered, so that ITPs can be nat- 
urally defined on macroscopic scales that remain small in 
comparison with the system size [ij. Beyond local equi- 
librium, more phenomenological endeavors based on ther- 
modynamical grounds ^ have been proposed, as well as 
statistical approaches illustrated on some specific mod- 
els [1, 0, [E B 13 ■ Finally, in the case of non-stationary 
slow dynamics, a notion of temperature may be derived 
from ageneralized fluctuation-dissipation relation (FDR) 
d, [§, llO, |Tl| in analogy to equilibrium statistical mechan- 
ics. Still, in spite of these numerous propositions, the 



relevance of this concept of effective temperature and its 
possible generality have been barely discussed. 

Besides, other notions of ITPs have appeared in the re- 
cent literature on nonequilibrium systems. For instance, 
in the context of stochastic models with a conserved mass 
M (or number of particles), like the Zero Range Process 
or different kinds of mass transport models jl2l |. a for- 
mal grand-canonical ensemble has been defined, in which 
systems with a total mass M appear with a probability 
weight proportional to exp(— /iM), /i being called a chem- 
ical potential. We call this ensemble a "formal" one, since 
no definition of the chemical potential is given prior to 
the grand-canonical construction, and no physical mech- 
anism allowing for fluctuations of the total mass (like, 
e.g., a contact with a reservoir) is described. Hence, 
the grand-canonical distributions introduced so-far ap- 
pear more as mathematical tools with interesting prop- 
erties, as it may be considered as a Laplace transform 
with respect to M of the canonical distribution for which 
M is fixed [l3| . Even more importantly, if this formal 
grand-canonical ensemble was to be considered as defin- 
ing the chemical potential, nothing could be said from 
it on the possible equalization of this parameter between 
two subsystems of a globally isolated system, since in 
the grand-canonical ensemble, the chemical potential is 
externally imposed. 

The aim of the present work is to introduce a pre- 
cise and general theoretical background allowing for the 
definition of ITPs conjugated to conserved quantities in 
nonequilibrium systems. Note that the systems consid- 
ered here are out of equilibrium not due to the presence 
of gradients imposed, for instance, by boundary reser- 
voirs, but because of the breaking of microreversibility 
(that is, time-reversal invariance) at the level of the mi- 
croscopic dynamics in the bulk. Accordingly, the ITPs 
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are not space dependent, as would be the case for systems 
that fulfill the local equilibrium assumption. We state 
the hypotheses underlying the present construction, and 
clarify the physical interpretation of the grand-canonical 
ensemble. We then discuss the relevance and usefulness 
of the concept of ITP, with particular emphasis on the 
description of phase coexistence. Although the proposed 
generalization of ITPs appears to be rather natural, it 
turns out that non-trivial problems arise as soon as sys- 
tems with different dynamics are put into contact. This 
issue is essential if one wants to measure an ITP using 
a small auxiliary system -just like temperature is mea- 
sured with a thermometer- and points at one of the main 
difficulties of nonequilibrium statistical mechanics. Note 
that part of the results reported here have appeared in a 
short version \XM- 

The paper is organized as follows. Sec. HIl introduces 
the concept of intensive thermodynamic parameters in 
the framework of nonequilibrium systems. The main 
condition of validity for this concept is given and some 
features in analogy to equilibrium statistical physics are 
discussed. Further the applicability of the definition for 
systems described with matrices (either matrix product 
ansatz or transfer matrix method) is studied, and an 
illustration of this approach on two simple models is 
given. Sec. [Till is dedicated to the issue of phase sepa- 
ration, specifically the problem of condensation, in dif- 
ferent models. We show how the concept of ITPs could 
improve the understanding of phase separation in these 
kind of models. Finally, the contact of two systems with 
different dynamics is investigated in Sec. lIVl with empha- 
sis on the issue of equalization of ITPs. We also explore 
possible ways of determining the ITPs in experiments or 
numerical simulations. 



II. NONEQUILIBRIUM INTENSIVE 
THERMODYNAMIC PARAMETERS 

A. Framework and definitions 

Let us start by considering a general macroscopic sys- 
tem that exhibits a steady state, and such that the dy- 
namics conserves some additive quantities, referred to as 
Qk, k = !,...,£, in the following. Such systems have 
been extensively studied for instance in the context of 
markovian stochastic models, and simple examples in- 
clude the zero-range process (ZRP) more general 
mass transport models [l^, [la, [13]; ot the asymmetric 
exclusion process (ASEP) in a closed geometry [l^. On 
the microscopic level, the nonequilibrium character of the 
dynamics manifests itself (apart from the lack of detailed 
balance) in the presence of a non-trivial dynamical weight 
fa associated with each microscopic configuration a, in 
the steady-state probability Pa- More precisely, the lat- 



ter reads 

^°- z(J:.q,) n/<«-Q>'- <" 

where the product of delta distributions ensures the con- 
servation of the quantities Qk- The function 

e 

Z{Qi,...,Q,)^J2f'=']l^^Qk-Qk) (2) 

a k=l 

serves as normalization factor, which will be referred to 
as "partition function" in analogy to equilibrium statis- 
tical mechanics. Let us emphasize that in equilibrium 
systems, the probability weight is either a constant 
independent of a, or an exponential factor accounting 
for the exchange of a conserved quantity with a reser- 
voir. For instance, in the equilibrium canonical ensemble, 
the conserved quantity would be the number of particles, 
whereas would be the Gibbs factor exp{—f3Ea), where 
Ea is the energy of configuration a, and (3 is the inverse 
temperature. On the contrary, in a nonequilibrium sys- 
tem, the weights /„ also account for purely dynamical 
effects related to the absence of microreversibility (the 
latter being deeply rooted in the hamiltonian properties 
of equilibrium systems), so that these weights generically 
depart from a constant, even if no conserved quantity is 
exchanged with a reservoir. 

To introduce a definition for ITPs in nonequilibrium 
situations, we first recall that their equilibrium defini- 
tion is related to the exchanges of conserved quantities 
between subsystems, and ensures the equality of ITPs in 
different parts of the system. Following the same line 
of thought in a nonequilibrium context, let us divide 
our system, in an arbitrary way, into two subsystems 
Sa and Sb- The sum Qka+Qkb — Qk is kept constant 
due to the conservation law, whereas exchanges of these 
quantities between the two subsystems are allowed. The 
microstate a is now defined as the combination of the 
two microstates {aa,Oib\ of the subsystems, so that the 
probability of a microstate is denoted as Pa^^at- ™" 
portant quantity in the following approach is the con- 
ditional probability \E'(Qia, • ■ • , Qea\Qi, ■ ■ ■ , Qi) that the 
conserved quantities have values Qka in subsystem 5q, 
given their total values Qk 

I 

*({Qfea}|{Q4) = E ^""."^ n ^^Qt: - Qka) ■ (3) 

OayOtb k=l 

The key assumption in the following derivation is that 
the logarithm of '^{{Qka}\{Qk}) satisfies an asymptotic 
additivity property, namely 

ln*({Qfca}|{Qfe}) = IriZaiiQka}) + lnZi,i{Qk - Qka}) 
- lnZ({Qfe}) + ejv({Qfca},{Qfc}) (4) 

with 

\eN{{Qka},{Qk})\ ^\\n^{{Qka}\{Qk})\ , (5) 
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in the thermodynamic limit N ^ oo. In Eq. ([4]), 
Zv{{Qku}) refers to the isolated subsystem S^,iyG {a, b}, 
and N is the number of degrees of freedom. That this ad- 
ditivity condition is fulfilled for some rather large classes 
of nonequilibrium systems will be illustrated in the fol- 
lowing examples. 

Some of the simplest systems that fulfill Eqs. ^ 
are lattice models with a factorized steady state distri- 
bution (where sites are labelled by i — 1,...,N, and 
a = {aj) 



Pa 



1 



N 



Z{Qi 



,Qi) 



fe=i 




Qk 



(6) 

in which case the term vanishes. Well-known exam- 
ples of models with factorized steady-states are for in- 
stance the ZRP [12J and other general mass transport 
models [l^, [l^ . Accordingly, the physical interpretation 
of the additivity condition given in Eqs. (H]), ([5]) is that, 
on large scale, the system behaves essentially as if the 
probability weight was factorized, although the genuine 
probability weight may not be factorized. 

As mentioned above, our aim is to define a parameter 
that takes equal values within the two (arbitrary) subsys- 
tems Sa and Sb- Guided by the equilibrium procedure, 
we consider the most probable value of Qka, denoted by 
Qka' which maximizes the probability ^'({Qfca}|{QA;})- 
This most probable value satisfies 



51n*({Qfc4|{Qfc}) 



dQ 



ka 



QL 



0. 



(7) 



Using the additivity condition (j4]), ([5]), we obtain the 
following relation 



d\nZa 



dQka 



d\nZb 



dQ 



kb 



(8) 



Hence, it is natural to define the ITP conjugated to 
the conserved quantity Q/c in a nonequilibrium system as 



Afc = 



dlnZ 
dQk 



fc = 1, . 



(9) 



Once the steady state is reached, this quantity equalizes 
in the two subsystems iSo and Sb due to Eq. and 
thus satisfies the basic requirement for the definition of 
an ITP. 

Actually, for the approach to be fully consistent, one 
has to check that the value of does not depend on 
the choice of the partition, as long as both subsystems 
remain macroscopic. We show in Appendix [X] that this 
is indeed the case, at least under the assumption (which 
is consistent with the additivity condition) that In Z is 
extensive. To achieve this result one shows that the ITP 
Afc obtained from the subsystems is equal to the global 
ITP defined on the whole system from Eq. ([9]), indepen- 
dently of the partition chosen (see Appendix . Thus 
from now on, we compute the ITP Afc using Eq. ^ for 
the whole system. 



B. Nonequilibrium "grand-canonical" ensemble 

Now that we have defined the notion of ITP, it is nat- 
ural to try to introduce a grand-canonical ensemble de- 
scribing a system in contact with a reservoir of conserved 
quantities, that imposes its values of {A^} to the system. 
Let us consider a partition of a large isolated system into 
two subsystems Sa and Sb, such that one of the macro- 
scopic systems is much smaller than the other. The larger 
subsystem, say Sb, serves as a reservoir for the {Qk}, 
leading to a natural definition of a "grand-canonical" 
distribution. This distribution is obtained by integrat- 
ing the "canonical" distribution ([T|) over the degrees of 
freedom of the reservoir 



-Pa. = 



fa 



Z{Qi,...,Qi) 



at, fc=l 



/a. 



ZiQi,...,Qi) 



Zb{{Qk-Qta}) 



(10) 



Note that in the above equation, we assume the factoriza- 
tion property faa,at = foafat, that is we neglect possi- 
ble boundary contributions, similarly to the assumptions 
often made in the equilibrium context. Expanding the 
logarithm of the partition function Zb for ^ Qk, 
one finds to leading order in Q^" 



\nZbi{Qk-Qt:}) ~ \nZbaQk})-Y. 



dlnZb 



fc=i 



dQ 



kb 



^ ka 



lnZ,({Q,.})-^AfeQ^," , (11) 



where we approximated 



dlnZb 



dQ 



kb 



dlnZb 



Qk 



dQ 



kb 



Ql 



fe=l 



A. 



(12) 



This is justified in the limit Na/N 0. Therefore we ob- 
tain for the "grand-canonical" distribution the following 
expression 



Pa. = ^ 



fa. 



exp(^-^Afeg-) 



(13) 



^(Ai,...,A,) 

which defines the grand-canonical partition function 
Z(Ai,...,A,) = ^/„„ exp(-^A,OL») ■ (14) 

aa \ k=l ) 

Interestingly, the cumulants < QkiQk-i ■ ■ ■ Qk„ >c may be 
expressed as a derivative of the logarithm of the grand- 
canonical partition function Z({Xk}), which thus appears 
as the associated generating function, as follows: 



<Qk,Qk,---Qk^ >c = (-1)" 

d d 



d 



dXki dXk2 dXk 



(15) 

■\nZ{{Xk}) 
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Note that some of the indices among fci , . . . , fc„ may be 
the same. This result generahzes the corresponding well- 
known equilibrium result [l^ . Accordingly, checking the 
validity of such relations in a given system does not 
show that microstates compatible with the constraints 
are equiprobable, contrary to what is sometimes im- 
plicitely assumed (e.g., in the context of granular matter 

til)- 



C. Applicability for nonfactorized steady states 

As stated above every nonequilibrium system that ex- 
hibits a product measure automatically fulfills the addi- 
tivity property given in Eqs. ([4]), ([5]), since cn vanishes 
in this case. The aim of this section is to discuss some 
typical cases for which the steady-state distribution does 
not factorize, but still satisfies the additivity condition. 
As a result, the ITP framework is relevant for such sys- 
tems. There are two important classes of models with 
nonfactorized steady states: systems where the station- 
ary state can be expressed by a matrix product ansatz, 
and models described by a transfer matrix. 



product ansatz: 



N / N 

^^^'^'^^ ^ Z{Q) ^'J\ ^^^'^'^ ^ ( ^ '^^ " ^ 



(18) 



where M{q) is a square matrix. Let us also introduce the 
matrix R{Q) through 



N / N 

R{Q) = / X{[dq^ M{q{)] 5\y^<l^^Q 
•' i=l \i=l 



(19) 



so that the normalizing factor Z{Q) = TrR{Q). 
This leads for the conditional probability distribution 
^iQalQ) to: 

^{Qa\Q)^^TT[RaiQa)Rb{Q-Qa)] (20) 

Loosely speaking, the additivity condition holds if the 
last factor Tr[Ra{Qa)RbiQ — Qa)] behaves essentially as 
Za{Qa)Zb{Q — Qa)- It IS shown in Appendix IB] that the 
additivity property ([4]), ([5]) is generically fulfilled for a 
system described by a matrix product ansatz with finite 
matrices. Whether it also holds for some classes of infi- 
nite matrices M{q) remains an open issue. 



1. Matrix product ansatz 



2. Models with, transfer matrices 



Practically speaking, considering a matrix product 
ansatz means that, in a one-dimensional system with pe- 
riodic boundary conditions, one expresses the probability 
weight fa in the form 



N 



/„ = Tr[]M„ 



(16) 



i=l 



where Mq,. is a (possibly infinite) matrix, Ui is the state 
of site i, and Tr is the trace operation over the matrices. 
For nonperiodic boundary conditions, a slightly different 
ansatz is used, namely 



N 



(17) 



where the vectors {W\ and \V) are determined by the 
boundary conditions (reflecting boundaries or injection of 
particles from a reservoir for instance). Matrix product 
ansatz have proved particularly useful in the context of 
the ASEP, where particles obeying an exclusion principle 
(that is, at most one particle per site is allowed) perform 
a biased stochastic motion on a one-dimensional lattice 

To test the additivity condition ^ for such mod- 
els, let us consider a generic lattice model with peri- 
odic boundary conditions. A variable qi is defined on 
each site i = 1, . . . iV, and we assume that the quantity 
Q = X^iLi 1i is conserved by the dynamics. The steady- 
state distribution is assumed to be described by a matrix 



Another matrix method that has become very popu- 
lar, in equilibrium as well as in nonequilibrium statistical 
physics, is the transfer matrix one. The main idea of this 
method is to formulate the partition function in terms of 
a product of a matrix, the so-called transfer matrix. 

To give an example of the application of this method 
for nonfactorized steady states let us consider, as previ- 
ously, a one-dimensional transport model on a ring, with 
a local variable qi on each site i, and such that the sum 
Q = X^iLi 1i is conserved. The variables qi may either be 
discrete or continuous. Let us now assume a steady-state 
distribution of the form 



P{{1^}) 



N / N \ 

ZlQ)\b{q^.q^+l)5\y^^q^~Q\ (21) 



with qN+i = qi, and g{qi,qi+i) is a symmetric function. 
The partition function Z{Q) is given by 

ZiQ) = [ ll[dq^ 9iq^,q^+l)] s(J2q^-Q) , (22) 

Note that a model of this type has been studied in [13], 
showing interesting nonequilibrium condensation proper- 
ties. Let us introduce the quantity SNiQ,qi,qN+i) de- 
fined as 

SNiQ,qi,qN+i) = J dq2...dqN (23) 

N / ^ \ 

i=l \i=l / 
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where qi is no longer identified with qn+i, contrary to 
Eq. dUl)- Then one has 



Z{Q)= / dqi SN{Q,qi,qi 
Jo 



(24) 



(25) 



and the distribution ip{Qa\Q) be written as 
^{Qa\Q) = dqi dqN^ + 1 

^ SNa{Qa,(ii,qNa+i) SnAQ ^ Qa,qNa+i,qi) 

^From this expression, a calculation similar in spirit to 
the one presented in Appendix [B] for the case of matrix 
product ansatz, allows one to show that the additivity 
condition ([5]) holds. The derivation makes use of 

the Laplace transform Sn{sj qi, qN+i) of Sn{Q, 91, Qn+i) 
with respect to Q, which can be written as a matrix prod- 
uct (hence the name transfer matrix method): 

Sn{s, qi,qN+l) - e--(9i-9«+i)/2 Tf (gi, g^+i) (26) 

where the transfer matrix Ts is defined by 

Ts{q,q')^9{q,q')e-^^'^+''">/' (27) 

Let us mention, here again, that the derivation of the ad- 
ditivity condition relies on some properties of the transfer 
matrix that are well established for finite matrices, but 
that might not be fulfilled in some cases for infinite ma- 
trices. 



To calculate the ITP corresponding to the conserved 
mass in the system we need to find the dependence of the 
partition function Z on M: 

Z{M) = / Y[[dm, mf-^] ,5 [ ^ m, - M j , (30) 

where the integrals are over the positive real axis. A sim- 
ple rescaling = XiM reveals the searched dependence: 

. N / N \ 

ziM) = AfST^i"'-! / Y[[dx,xr']s[j2^^ - 1 

i=l / 



(31) 



with 1] = N ^ Vij ^i^d where Kpf is a constant in- 

dependent of M. The ITP is obtained from the derivative 
of InZ 



A 



d\nZ Nt]~1 



dM M 
leading in the thermodynamic limit — > oo to 



where p ~ M/N denotes the average density. 



(32) 



(33) 



D. ITPs at work on simple models 

1. Mass transport model with factorized steady-state 

As a first example we consider a simple one- 
dimensional mass transport model on a ring with only 
one globally conserved quantity, referred to as mass 
M = X^iLi ^ introduced in Ref. [Tsj . The masses rrii 
are a priori positive and real variables. The continuous 
time stochastic dynamics is defined as follows. A mass fi 
is transferred from a randomly chosen site i, containing 
the mass rrii, to site i + 1 according to the following rate 

^,i^,\m,)=vi^,)^^^^P^. (28) 

Thus transport is totally biased, which generates a flux 
of mass along the ring. With the above rate, the steady- 
state distribution is of the form [l^ (see also Appendix [Cl 
for a more general case) 

N / N \ 

Pi{rn^})^^^^^Mm.)5\Y2m.~M\ , (29) 

where the single-site weight fi{mi) may be site- 
dependent. Let us here consider for fi{mi) the simple 
form fi{mi) = rn^'~ with ?7i > for all i. 



2. Model with pair-factorized steady-state 

Let us now consider a second example to illustrate that 
our approach works as well on systems that do not exhibit 
a factorized steady state. We therefore consider a model 
similar to that used in the preceding section, but with 
a transport rate that depends, besides the mass on the 
concerned site, also on the masses on the two neighboring 
sites: 

/ I N / g{mi-i,m., - p) 

ip(p\m,^i,mi,mi+i) = v(p) ^ — x 

g(mj_i,TOi) 

g{m, - p.,m,+i) ^^^^ 

g{mi,mi+i) 

Note that this is a generalization of the dynamics in- 
troduced in defined for discrete masses. It can be 
shown (see Appendix [C| that these dynamics lead to a 
pair factorized steady state of the form 

Choosing, as a simple example, the function g{m, n) as 
g{m, n) = {m^n^ -f m.^n°'y (36) 
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with a, /3, 7 > 0, we obtain the behavior of the partition 
function on M, using again a simple scahng argument: 



Z{M) 



N / ^ \ 

Y\_[dmig{mi,mi+i)]5 - M 

1=1 \i=l / 

N 



{xi, a^i+i)] 



(37) 



with a prefactor K]\j independent of M. Thus the ITP 
conjugated to the conserved mass reads 

dlnZ jV[7(a + /3) + l]-l 

^ ~ ^ M 

yielding in the thermodynamic limit 

7(a + /3) + 1 



A 



(39) 



Note that here, we did not checked the additivity con- 
dition ([5]), while the distribution ((35|) would lead to 
an infinite transfer matrix. Yet, we check in Sec. IIV lT2l 
on a specific example, that the above results are indeed 
consistent. 



E. Discussion on the physical interpretation of 
ITPs 



Now that the formalism is introduced, let us insist on 
the physical interpretation of ITPs. First, it is impor- 
tant to note that, in the present framework, the grand- 
canonical ensemble is explicitely derived from the canoni- 
cal one, and thereby gains a clear physical status, whereas 
the "formal" grand-canonical ensemble often considered 
in the literature, is in a sense more like a Laplace trans- 
form, without physical ground . A nice illustration of 
this difference appears when considering the topology of 
the physical system which may indeed differ in the canon- 
ical and grand-canonical ensemble. For instance, when 
dealing with a canonical system on a one-dimensional 
ring, the corresponding grand-canonical system would be 
on a segment, the complementary segment playing the 
role of the reservoir. Such a subtility does not emerge 
when dealing with the Laplace transform of the canoni- 
cal distribution. Second, in the present context the ITP 
has the status of a well-defined thermodynamical vari- 
able, which equalizes between subsystems. We shall il- 
lustrate in the next section how it applies to the descrip- 
tion of phase coexistence and how it allows to extend 
our "thermodynamical way of thinking" to some out-of- 
equilibrium situations. 

One might also argue that, at least when the conserved 
quantity is a number of particles, the present defini- 
tion of the associated ITP may be recovered by mapping 



the nonequilibrium stationary distribution on an effective 
equilibrium canonical distribution (thus introducing an 
effective hamiltonian) and computing the corresponding 
equilibrium chemical potential. Although it leads to the 
same result as ours, such a procedure is less general and 
rather confusing, since the nonequilibrium aspects of the 
system are somehow hidden. Indeed, one may conclude 
from such a mapping that two nonequilibrium systems 
in contact should equilibrate their chemical potentials, 
which is not necessarily true, as will appear more clearly 
in Sec. [iVl 

Besides, it is sometimes thought that ITPs somehow 
identify with fluxes within a given system, for instance 
the fiux of mass between two neighboring sites in a mass 
transport model. One reason for this is that in one- 
dimensional models, the stationary flux necessarily takes 
the same value throughout the system. Yet, this identi- 
fication is actually not valid in general, but only for spe- 
cial cases such as for instance the one-dimensional ZRP, 
where the flux $ of particles within the systems is equal 
to the fugacity e~^. To illustrate this point, let us take 
for example a simple mass transport model on a ring with 
N sites and transport rates (p{fi\in) defined in Eq. (|28p . 
assuming for simplicity that fi{m) = f(jn) for all i. The 
average fiux (f> crossing a given link reads 



dmP{m) I dfj. fi(p{fi\m) 

^0 

Using the single site probability P{m) — Cf{m)e~^ 
being a normalization constant, one finds 



(40) 
C 



= C/ d^i^iv{^l)e-^^' dxf{x)e-^'' (41) 

Jo JO 

with the change of variables x = m ^ fi. Using the defi- 
nition of C, the above equation reduces to 



(42) 



Considering this equation we can conclude that the rela- 
tion between the flux and the ITP can be highly nontriv- 
ial, and that the two notions should not be identifled. 
Moreover, there is no obvious reason why, in a generic 
system of dimension d > 1, the local flux should take 
equal values throughout the system. Note also that in 
models with closed boundary conditions, the global flux 
of a conserved quantity is forced to be zero, whereas ITPs 
a priori take nonzero values, making them a more useful 
characterization of non-equilibrium steady states. 



III. RELEVANCE OF ITP'S TO DESCRIBE 
PHASE COEXISTENCE 

A. Principle of the approach 

One of the main interest of the notion of ITP, which 
made its success in equilibrium, is that such parame- 
ters take equal values in different subsystems of a given 
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system, regardless of their macroscopic state. Indeed, 
subsystems may for instance have different densities of a 
given conserved quantity, but they should have the same 
value of the associated ITP. A case of great interest where 
this situation arises is that of phase coexistence. At equi- 
librium, phase coexistence is described by the equality 
of the different ITPs (temperature, pressure, chemical 
potential). We shall now argue that the present ITP 
formalism allows for a similar description of nonequilib- 
rium phase coexistence, at least when this phenomenon 
is related to a conservation law, and when the additivity 
condition (g]), (O holds. 

As a simple illustration, we shall consider in this sec- 
tion a well-studied example of nonequilibrium phase co- 
existence, namely condensation transitions. Such transi- 
tions have been reported in ZRP |,12. 22]. [gsl. [23 | , and in 
more general mass transport models [la . [l7| . These mod- 
els have in common a critical density p'^"^ above which a 
condensation transition occurs, that is, a finite fraction 
of the total mass condenses onto a given site (or on a 
small domain 

The standard way to compute the critical density in 
ZRP for instance, is to use the grand-canonical partition 
function Z, considered as a function of the fugacity z = 
e-^ [11 : 



Z{z) = ^""ZMiM) 



N=0 



where z is fixed by the density 

_ M _ z d\nZ 
iV ^ N dz 



(43) 



(44) 



and to look for the convergence radius of Z{z) in the 
complex plane of z. To study the condensed phase in 
more details, it is necessary to use a canonical ensemble 
approach where the total mass is fixed [l^ . Accordingly, 
the present ITP formalism turns out to be well-suited for 
such a study. 

In Sec. IIIIBI and IIII CI we illustrate in a pedagogical 
manner how ITPs may give a natural quantitative de- 
scription of the condensation phenomenon, on the exam- 
ple of simple mass transport models. Qualitatively, the 
general procedure proposed is the following. Interpreting 
the condensation as the coexistence of a fluid phase and a 
condensed phase, one concludes that from the definition 
of ITPs, the value of the ITPs should equalize in the two 
phases. The ITP for the single-site condensate is often 
easily obtained, in which case the value of the ITP in the 
fluid phase is also known. Then the equation of state of 
the fluid phase, computed in the grand-canonical ensem- 
ble, can be used to determine the density of this phase 
(which turns out to be the critical density). Hence, the 
total mass of the fluid phase is known. ^From the knowl- 
edge of the total mass, one finally deduces the mass of 
the condensate. 

Such a description of phase coexistence is a good il- 
lustration of the application of ITPs to a rather simple 



out-of-equilibrium situation. It is also a first step to- 
ward the description of the more complex situation of 
the contact between two different systems. As we shall 
see in Sec. IIV[ this case may reveal some difficulties, un- 
expected within a formal analogy with equilibrium (for 
instance by defining a formal grand-canonical distribu- 
tion, or through the introduction of an effective hamilto- 
nian), but well enlightened in the present framework. 



B. Mass transport model with one species 

1. Homogeneous model 

As a first example of the application of the ITP concept 
in phase coexistence let us return to the model described 
in chapter III D 11 defined by the transport rates given in 
Eq. (|28p . Assuming that fi{m) = f{m) for all z, and that 



/(to) 



for TO ^ cx), one finds for 7 > 2 a condensa- 



tion if the average density exceeds a critical value ||12i.il6l|. 
In the following we revisit this condensation transition, 
and show how our approach of equalized ITPs can reveal 
the physics behind this phenomenon. If we assume that 
the condensate, which occurs on a randomly chosen site 
J, carries the macroscopic mass Mc = 0{N), its canoni- 
cal partition function yields 

Zc{Mc) = / dnij f{mj)S{mj - Mc) 
Jo 

= f{Mc) ^ M--* . (45) 

The ITP conjugated to Mc is the chemical potential Ac 
of the condensate, which is given by 



Ac 



din Zc 

dMr 



7 
' Mc 



(46) 



In the thermodynamic limit Mc — ^ 00, one thus obtains 
Ac = 0. The equality of ITPs for the condensed and 
the fluid phase therefore leads to A/ = Ac = 0. The 
density p/ of the fluid phase, associated to A/, can be 
determined from the equation of state computed in the 
grand-canonical ensemble 



Pf 



dm mf (111)6 ^f"^ dmmf{m) 
dm f{m)e^'^f™ dm f{m) 



(47) 



and exactly gives the critical density p'^"* . Note that 
for a value of 7 < 2 the critical density becomes infinite, 
which means that no condensation occurs. Accordingly, 
the mass of the condensate is given by 



Mc = M - Np" 



(48) 



if the overall density p — M/N is larger than p'^"*. This 
therefore leads to a thorough description of the conden- 
sation in this system. 
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2. Model with an impurity 



C. Mass transport model with two species 



Another well-known situation where condensation oc- 
curs is when a single impurity site exhibits a dynamics 
which differs from those of the other sites [l^. In this 
case the condensation no longer occurs on a randomly 
chosen site, but on the impurity itself. One of the sim- 
plest choice for the weights in such a model corresponds 
to 



fi{m) = /imp(m) 



with C > 



(49) 



for the impurity site and 

Mm) = /hom(m) = m'!-^ with ?7 > (50) 

for the remaining sites i > 1 ('hom' stands for 'homoge- 
neous'). In this case the canonical partition function for 
the condensate reads 



which amounts to an ITP for the condensate 

din Zr 



A. 



dMr 



(51) 



(52) 



This yields for the fluid phase A/ = ^, again by equalizing 
the ITPs of the two phases. The equation of state Pf{Xf) 
for the fluid phase, computed in the grand-canonical en- 
semble, reads 



/p°°dTOm''-ie-^/™ ~ A/r(7y) ^ A/ 



(53) 



Thus the equality of the chemical potentials forces the 
fluid to have a fixed density po — ^/d long as the 
condensate is present, that is for p — M/N > pQ. Inter- 
estingly, po is not in itself the maximum density of the 
homogeneous fluid phase, but is simply a density imposed 
by the impurity. 

Finally, let us also briefly mention the interesting case 
where the impurity is defined by 



/imp(m-) 



(54) 



while the homogeneous phase is still defined by Eq. ([SC 
Assuming again that there is a condensate on the im- 
purity, one finds Ac — 2C,Mc — + oo. Thus one expects 
A/ — > oo, yielding Pf — Q using Eq. ([53]) . This case is of 
particular interest, since on the one hand there is a full 
localization of mass on the impurity in the thermody- 
namic limit, and on the other hand it cannot be studied 
within the framework of the grand-canonical ensemble, 
due to the faster-than-exponential divergence of /i(r7i). 

Accordingly, the above calculations provide a simple 
description of the condensation in terms of the (out-of- 
equilibrium) "equilibration" of two coexisting phases. In 
particular, that the condensate may absorb all the excess 
mass is understood due to the fact that Ac is independent 
of Mc. 



The approach explained above for a model with one 
conserved quantity can easily be extended to more com- 
plicated dynamics where more than one species are in- 
volved. In the case of two species, the general strategy is 
the following. Let us consider a lattice model with two 

conserved masses Mi = Yl^=i and M2 = J2iLi "^2i- 
We assume that a condensate takes place on a given site 
io. The key point of the approach is that the ITPs con- 
jugated to Ml and M2 take equal values in both phases. 
We shall denote these common values as A^, k = 1,2, 
regardless of the phase we are dealing with. Assuming a 
factorized steady-state, one can easily compute the values 
of Afc in the condensate as a function of the condensed 
masses Mic and M2C, leading to two equations of the 
form: 



Afe = Afec(Mic,M2c) , 



1,2 



(55) 



By analogy with the single-species case, it is expected 
that at least one of the ITPs vanishes in the condensate in 
the thermodynamic limit, as long as we consider a system 
without impurities; we thus assume that Ai = 0. Then, 
one needs to compute the densities in the fluid phase us- 
ing the grand-canonical equation of state pk / (Ai = 0, A2). 
Finally, the conservation of the total mass generically 
yields two coupled non-linear equations for the unknown 
variables Mic and M2c: 

Mfec + NpkfiO, A2c(Mic, M2c)) = Mfe , A: = 1, 2 . (56) 

These equations may be hard to solve analytically in gen- 
eral, but they could be solved numerically. In addition, 
simplifications may appear in some cases, as we now il- 
lustrate on a specific model. 

Let us consider a model with two conserved masses 
as above, which are being transported on a lattice of 
arbitrary dimension d, with periodic boundary condi- 
tions. The transfer rates for an amount of mass fii or 
P2 from a random site to a neighboring one are given by 
ipi{p.i \mii,m2i) or </?2(M2|™iij ra2i) respectively, in accor- 
dance with 



ipi{pLi\mii,m2i) 

<P2(M2|"ll^,"^2^) 



h{mu - p.i,m2i) , . 
/i(mii,TO2i) 

V2[P2) r? ^ — (58) 



h{mu,m2i) 



(59) 



For lattices of dimensions d > 1, one allows a symmetric 
mass transfer in any direction transverse to the flux with 
the same rates ^i{pi\'mii,'m2i) and ^P2{p2\'mii,m2i) as in 
the direction of the flux. The steady-state distribution 
in arbitrary dimension d is then given by 



Z{Mi,M2) 



ruki — Mk 



k=l 



(60) 
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Guided by the studies of the two-species ZRP [2j 
choose an exponential form for the weights 



/i(mi, 7712) = exp(- 



Kmim2 



(61) 



with K, (T > 0. To obtain the conditions under which a 
condensation occurs, we calculate the canonical partition 
function of the single-site condensate Zc{Mic, M2c) = 
h{Mic, M2c)- Hence one obtains for the chemical poten- 
tials of the two species, computed in the condensate 



Ai = 



d In Zc 
dM2c 



M. 



K<7 



2c 

Ml, 



(62) 
(63) 



Therefore in the thermodynamic limit, one finds for the 
chemical potential of the first species Ai = 0. The den- 
sities in the fluid phase arc then computed in the grand- 
canonical ensemble approach, which yields for pif: 



Pi/(0, A2) 



dm2 6-^2^2 dmi e-™i™2 " 
1 r{l + 2a) 



(64) 



For p2 f , one obtains in the same way 

Jg°° dm2 m2e^^'^"^'^ drrii e~ 



P2/(0, A2) 



(65) 

The integration over mi is straightforward and leads to 

1 r(2 + cr) (T + 1 



P2f 



\2T{1 



A. 



(66) 



Hence if a condensate exists, one has A2 > to en- 
sure a finite positive density p2f in Eq. (jM]) . Then, 

from Eq. §^ one finds that M2C Mi/^''+^\ which 
means that the condensed mass M2C behaves subexten- 
sively, and the density of the fluid phase for the second 
species is thus equal to the overall density, p2f = P2- 
Due to Eq. ([66]) it follows that A2 is determined only by 
the value of p2- Therefore, the value of pif is also fixed 
by the overall density p2 of the second species, since it 
is only dependent on A2, as seen in Eq. As a re- 

sult, a condensate forms only under the condition that 
Pi > pI"\p2) given by 



pr*(P2)=Pi/(0,A2(p2)) 



r(l + 2cr) / a + 1 
kT{1 + (7) \ P2 



(67) 

Using Eqs. ([55]) and as well as the relation Mic — 

N{pi — **), one finds the following expression for the 
condensed mass of the second species 



M2e = 



K<7p2 



Nip, - pr') 



(68) 



To sum up, we have shown in this section that in the 
framework of ITPs, a simple procedure can be developed 
to describe the coexistence of different nonequilibrium 
phases in mass transport models, based on the idea that 
ITPs equalize their values in these phases. We now turn 
to the more complex situation, where systems with a dif- 
ferent microscopic dynamics are put into contact. 



IV. CONTACT OF SYSTEMS WITH 
DIFFERENT DYNAMICS AND ITP MEASURE 

A. General approach 

Within the equilibrium context, the temperature, pres- 
sure and chemical potential of two different systems put 
into contact equalize, as long as the contact allows the 
conjugated conserved quantity (energy, volume or par- 
ticles) to be exchanged. This is actually a strong state- 
ment, as it is true even if the two systems considered have 
very different microscopic dynamics, provided that both 
systems can be described by an hamiltonian. Whether 
such property also holds in the present more general con- 
text of nonequilibrium ITP formalism is thus an essential 
issue. Potentially interesting applications are the descrip- 
tion of the effect of the environment on a system, or the 
possibility to measure the value of an ITP using a small 
auxiliary system, in the same way as temperature is mea- 
sured with a thermometer. 

As we shall see in this section, the equalization of the 
values of the ITPs in two systems in contact is actually 
not automatically fulfilled in a nonequilibrium context 
due to the necessity to satisfy the additivity condition 
dH), (O, and the way the two systems are connected has 
to be considered carefully. To highlight this point, let us 
examine two different systems (that is, with different mi- 
croscopic dynamics) iSi and ^2 that separately conserve 
the same physical quantities Q, with values Qi and Q2 
respectively. When put into contact, the dynamics at 
the interface generates the distribution ^{Qi\Q) for the 
random partition of Q = 9* i^^^o Qi (in system Si) 
and Q2 — Q ~ Qi (in system ^2) respectively. Assum- 
ing, in the spirit of equilibrium calculations, that the two 
systems are weakly coupled, i.e. the only coupling is the 
exchange of Q, the global probability distribution reads 

rQ 

Pi{q^}) = / dQi<^>{Qi\Q)Pi{{q,,}) P2i{q,,}) 



Zi{J2i,eS, 9n)^2(Q - EiiG5i In) 

Fl{{q^A)F2{{q^,})^[J2q^-Q] (69) 



where Fi{{qi^ }) and ^2 ({<Zi2 }) the probability weights 
of systems Si and 52 respectively, taken as isolated. For 
the ITP to equalize, the additivity condition (jH), ^ has 
to be satisfied, when applied to Si and 52 considered as 



10 



the two subsystems of the global system. Intuitively, this 
additivity condition means that the probability weight es- 
sentially factorizes, as already mentioned when describ- 
ing the grand-canonical ensemble. Due to the appear- 
ance of the quantity X^iieSi '^^i (ESJi nonlocal con- 
tributions arise in the probability weights appearing in 
P{{qi}). For these weights to factorize (and thus for the 
additivity condition to be fulfilled), the prefactor ^/ZiZ2 
has to be equal to a constant (since it would otherwise de- 
pend on the nonlocal quantity X)iie5i In) which is noth- 
ing but l/Z{Q), up to corrections vanishing in the ther- 
modynamic limit. Thus <&(Qi|Q) should necessarily be 
of the form 

(70) 

again up to possible small corrections. This result may 
alternatively be interpreted in the following way. The dis- 
tribution $(Qi IQ) is actually nothing but the conditional 
distribution '^{Qi\Q) introduced in Eq. for the spe- 
cific partition of the global system into (iSi, 1S2) -whereas 
^(QilQ) is a priori defined for an arbitrary partition. 
Eq. ([70| is then simply the additivity condition (jl]), ^ 
applied to the partition {Si,S2)- 

When putting into contact two nonequilibrium sys- 
tems, ^{Qi\Q) does not obey Eq. ([70)) in general, so that 
ITPs do not necessarily equalize. Yet for special cases, 
this equalization may be recovered, as we shall see in 
Sec. IIVBI In any case, it is an important challenge to be 
able to understand which microscopic properties the dy- 
namics of the contact has to satisfy so that equalization 
of ITPs hold. This is the topic of Sec. HVCl 

To clarify the relation with the equilibrium case, let us 
consider for instance the equilibrium canonical ensem- 
ble, where the conserved quantity would be the number 
of particles. The system is described by a Boltzmann- 
Gibbs distribution, and assuming that the hamiltonian 
does not include long-range interactions, the additivity 
condition (g]), ([5]) is satisfied. ^From this, one deduces 
that Eq. ([70l) is necessarily satisfied, since otherwise the 
additivity condition would not hold. 

B. Connecting two different mass transport models 

To illustrate an implementation of the above idea, we 
consider two single-species mass transport models in con- 
tact, as represented schematically in Fig.[TJ The dynami- 
cal rules are site-independent within each system, but are 
different in iSi and S2- To be more specific, the transport 
rate within system Si, reads 

^^{^\m) = v{^)Mp^ , ^=1,2. (71) 

Let us emphasize that the function v{fi) is the same for 
both systems. The reason for this choice, ensuring the 
equalization of the ITPs, will be explained in detail in 
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FIG. 1: Illustration of two mass transport models with peri- 
odic boundary conditions, connected through a contact made 
of a single additional link. A zoom over the contact is also 
shown in the dashed frame; arrows indicate the oriented links 
through which mass is transported. 

Sec. IIV CI A particular case where this condition holds 
is the ZRP one, where /i = 1 is the only allowed value, 
and v{l) is set to 1. 

In order to define the dynamics at the contact, we 
extend the dynamical rules in such a way that a mass 
located at site ii £ Si at the interface (see Fig. [T|), is 
transferred with rate ipi{/j,\m) either to the neighboring 
site + 1 or with the same rate to site 12 G ^2. A similar 
rule is applied for a mass located on site 12, regarding the 
transfer to zi. Hence the two homogeneous systems sim- 
ply combine to an inhomogeneous one, for which one can 
find a factorized steady-state, as v{^) is site-independent, 
and as the graph on which the system is defined satisfies 
the required geometrical constraint [IHl- It follows that 
in this case the additivity condition ([4]), ([5]) holds, and 
the ITPs of the two systems equalize. 

Let us now give some examples illustrating the conse- 
quences of the equalization of ITPs for two systems in 
contact, within the framework of the above model. It 
should first be noticed that the equalization of ITPs en- 
forces constraints on the densities of the two systems, 
which may thus differ one from the other if the two sys- 
tems have different equations of state. An interesting 
situation arises in the presence of a condensate. To be 
more specific, let us consider two different mass transport 
models, without impurities. The two systems are initially 
both isolated, and we assume that the first one contains 
a condensate, while the second one does not. This means 
that, in this initial stage, Ai = whereas A2 > 0. When 
put into contact, the dynamics will be such as to equalize 
the ITPs Ai and A2, through a transfer of mass between 
the two systems. As A generically decreases when the 
density increases, mass has to be transferred from Si to 
^2. Withdrawing mass from Si actually does not affect 
the fluid phase of Si (at least in a quasi-static limit) as 
long as the average density pi remains above the critical 
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FIG. 2: Time dependence of the masses of the condensates 
(a), the ITPs (b) and the densities of the fluid phases (c) in 
two mass transport models in contact, 5i (sohd line) and <S2 
(dashed line) -see Fig. [1] Both of them contain an impurity 
site. Parameter values are = 1.5 and ^2 = 1 for the im- 
purities, rji = 2 and 772 = 1 for the homogeneous parts -see 
text. The systems are separately in steady state (with overall 
initial densities pi — p2 ~ 2) before the contact is established 
at time t=0 (dotted line). Note that the ITPs eventually 
equalize, and that the difference between the fluid densities 
increases after the contact is established. 



density pi"* , and mass is taken from the condensate only, 
in this first stage. Assuming that the critical density of 
S2, pi*^**, is infinite, the condensate in iSi eventually dis- 
appears, and both ITPs converge to a strictly positive 
value that hes between the initial values of Ai and A2. 
Yet, the final densities of the two systems are different in 
general. 

One may also think of more complex situations. Let us 
consider the case where each of the two systems contains 
an impurity (as defined in Sec. IIII B 2|) . with respective 
parameters ^1 and ^2 (C2 < ■Ci)- addition, we con- 
sider different functions fhom{m) for the homogeneous 
part of each system, namely f\iom,u{m) = m''""^ in sys- 
tem iSjy, V = 1,2. Before the contact is switched on, 
both systems are in steady state and contain a conden- 
sate, so that Ai = ^1 and A2 = ^2- When the contact 
is established at t = 0, the ITPs tend to equalize, and 
as A2 < Ai, mass is transferred from 1S2 to iSi (as p de- 
creases with A), independently of the values of rji and 772. 
But as iSi already contains a condensate, the density of 
the fluid phase is fixed, and all the mass brought to Si is 
actually transferred to the condensate. Hence, the final 
state corresponds to Ai = A2 = ^i, with Si containing a 
larger condensate than initially, and ^2 being at a den- 
sity P2(A2 = 6) < P2'** = P2(A2 = 6), so that there 
is no more condensate in 1S2 . Fig. [2] presents the results 
of numerical simulations of the above situation, showing 



in particular that the condensate of ^2 disappears once 
the contact is established. In addition, this simulation 
confirms that the ITPs of the two systems are control- 
ling the direction of the flux of mass. Indeed, mass is 
transferred from 1S2 to Si even though the density of the 
fluid is larger in Si than in 52, which might seem rather 
counterintuitive given that the two systems are in contact 
through their fluid phases only. 

Note also that strictly speaking, the determination of A 
using the equation of state is valid only in a steady state, 
so that one should wait until the density p is stationary 
before determining A from p. Still, for the purpose of 
illustration, we present here \{t) deduced from the fluid 
density Pf{t) even in the nonstationary regime. This is 
meaningful if the exchanges between the two systems are 
slow enough so that each system may be considered in a 
quasi-steady state. 

^From the above result, one sees that the impurity in 
system iSi (that is, the one with the larger value of ^) 
plays the role of a reservoir of mass that fixes the value of 
the ITP of ^2 to A2 = ^2- Interestingly, Si only needs to 
contain a mass of the same order as that of 1S2 to act as a 
reservoir, as the ITP Aic of the condensate is independent 
of its mass Mic as long as Mic is macroscopic. On the 
contrary, usual reservoirs require a mass much larger than 
that of the systems they are in contact with. 

As we have seen on these simple examples, the notion 
of ITP allows one to make, essentially without calcula- 
tions, non-trivial predictions about the behavior of mass 
transport models put into contact (for instance, the fi- 
nal densities of both systems); only the knowledge of the 
equation of state for each system taken separately is re- 
quired. However, such predictions can be made only if 
the ITPs equalize. In the present case of mass transport 
models, this equalization holds as long as the transition 
rates obey the relation (ffTj) . What are the conditions for 
such an equalization to hold in more general situations 
will be the topic of the next section. 



C. Characterization of the dynamics at the contact 

In this section, we consider the more general case of two 
systems in contact as schematically illustrated in Fig. [3l 
The dimension d of the two systems is arbitrary, and the 
contact consists in a set of links between the two systems. 
Generically, the contact relates some parts of the borders 
of each system, as shown in Fig. [3l but one may also think 
of more complex types of contact. 

As mentioned in Sec. IIV AI the dynamics at the con- 
tact plays an essential role in the possibility to equalize 
the ITPs of two connected systems. This equalization oc- 
curs if the condition (|70p . or equivalently, the additivity 
condition ([5]), is fulfilled. Yet, testing this condition 
is a priori very difficult, and its interpretation in terms of 
the microscopic dynamics at the contact is not obvious. 
On the other hand, if two systems Si and ^2 are con- 
nected and globally isolated, the flux 0i_+2 transferred 
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FIG. 3: Schematic drawing of two generic systems in contact. 
Arrows illustrate the possibility to transfer a conserved quan- 
tity between the two systems. Contrary to Fig.[l] the contact 
generally extends over many sites. 

from Si to 1S2 through the contact has to be equal, in 
steady state, to the reverse flux (f>2^i going from S-? to 
Si. 

For the sake of simplicity, we now introduce some as- 
sumptions that allow one to analyze, at least in a simple 
case, this generic problem of the contact between two sys- 
tems. For definiteness, let us consider lattice systems Si 
and 52 that may exchange a globally conserved quantity. 
The set of sites belonging to the contact in system S,j, 
V = 1, 2, is denoted as C^. As already mentioned, the sit- 
uations considered correspond to the weak coupling limit. 
In the present context, this means that the flux per site 
crossing the contact is typically much smaller than the 
local flux between two sites of a given system, so that the 
contact does not perturb the dynamics of each system, 
apart from the (slow) exchange of Q. In addition, the 
specific assumptions used in the following arguments are 
that: 

(i) the flux 01^2 depends only on Ai, and not on the 
properties of 52 such as A2; respectively, (/)2^i depends 
only on A2 (yet, the total flux <^i^2 — 4>2^i depends on 
both Ai and A2); 

(ii) the probability weights of Si and 1S2 , each one con- 
sidered as isolated, are factorized as products of one-site 
weights. 

In the spirit of hypothesis (i), the dynamics at the 
contact is defined by the probability rate 93^^ (/ijgij ) to 
transfer ji from site ii G Ci in iSi to ^2; a similar rate 
^i^{^A1i2) defines the transfer from 12 £ C2 in ^2 to Si. 
Under these assumptions, one can compute the fluxes 
0i^2(Ai) and (/'2-*i(A2); in particular, 0i^2(Ai) reads 

roo pqt 

<^i^2(Ai) = ^/ dq,Pi{q,) d^l^Ji'p'i[^J.\qi). (72) 

where Pi{qi) is the single site probability distribution in 
Si. Given that Pi{q) = cifi{q)eicp{—Xiq), ci being a 
normalization constant, it follows that 

0i^2(Ai) = ci dq d^ifi^[°\fi\q)fiiq)e-^''i , 
Jo Jo 

(73) 

where we have introduced the "total" rate 

ieCi 



By exchanging the indexes of the systems, a similar rela- 
tion holds for (l)2^i{^2). 

For the two systems to equalize their ITPs, it is nec- 
essary that the equality 0i^2(Ai) = (/'2^i(A2) leads to 
Ai = A2, for arbitrary values of (say) Ai. In other words, 
the two functions (j)i^2i^) and 02^i(A) must be iden- 
tical. Let us then compute the difference 0i_>2(A) — 
02-^1 (A). After a straightforward calculation, this dif- 
ference can be expressed as 

/•oo /"OO /•qi 

4>l^2W - 4>2^lW = C1C2 / dqi / dq2 / dfi 

Jo Jo Jo 

X Aie-^(«^+«^)[<^r(Mki)/i(gi)/2(g2) 

-vTif^\Q2 + - M)/2(g2 + m)] (75) 

In order that (/)i^2(A) — 02^i(A) vanishes for any value 
of A, it is necessary and sufficient that the expression 
between brackets vanishes, that is 

<*(Mki)/i('?i)/2(g2) (76) 
= </'2°*(Ai|92 + m)/i(9i - m) 72(92 + fJ.) 

One can check that Eq. ([7S)l is precisely a detailed bal- 
ance relation between configurations (91,(72) and {q'l = 
qi — ^, ^2 = 92 + A*)- Yet, let us emphasize that this 
detailed balance relation does not concern the true mi- 
croscopic dynamics at the contact, but rather an effec- 
tive, coarse-grained dynamics, defined by (/3^°*(/i|9) that 
reduces the contact to a single effective link. To illustrate 
this point, one can imagine a contact made of two fully 
biased links between Si and 52. The first link can only 
transfer the conserved quantity Q from Si to ^2 , whereas 
the second one only allows Q to be transferred from ^2 
to Si. Then, the microscopic dynamics at the contact 
does not satisfy detailed balance, but the coarse-grained 
dynamics may fulfill this condition. 

Note that if the true contact already consists in a sin- 
gle link as illustrated on Fig. [U the effective dynamics is 
the true one, so that the dynamics at the contact really 
satisfies detailed balance. This is basically the interpre- 
tation of the condition used in Eq. (|7T|) . according to 
which w(/i) has to be the same in the two systems: this 
condition ensures detailed balance at the contact, even 
though detailed balance breaks down in each system due 
to the presence of fluxes. In constrast, if the functions 
f (/i) were to be different in the two systems the ITPs 
would not equalize due to the lack of detailed balance at 
the contact. 

In summary, we have seen the importance of the dy- 
namics of the contact, which must fulfill condition ([70)1 to 
ensure the equalization of the ITPs among the connected 
systems. On the other hand, the effective detailed bal- 
ance Eq. ([75)1 is a priori a less general statement, due to 
the simplifying assumptions used to derive it. However 
it provides us with a physical interpretation of the con- 
ditions required for the dynamics at the contact. One 
may hope that the resulting physical picture might be 
relevant beyond the strict validity of Eq. (|76p. 
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D. Measure of an ITP 

1. Notion of I TP-meter 

An essential issue about ITPs would be the ability to 
measure them. In analogy to the measurement of tem- 
perature in equilibrium statistical mechanics, one could 
imagine to realize such a measurement by connecting an 
auxiliary system to the system under consideration. This 
problem has also been addressed in In this general 
framework, we shall call this (perhaps conceptual) in- 
strument "ITP-meter" , in reference to the nomenclature 
"thermometer" . There are three essential requirements 
which have to be fulfilled to perform such an implemen- 
tation. First, the values of A in the gauged system and in 
the ITP-meter need to equalize. Second, the measured 
system must not be disturbed by the measurement, so 
that the ITP-meter should be small with respect to the 
system upon which the measure is performed. Third, 
the equation of state of the ITP-meter has to be known, 
since one needs to deduce the value of its ITP from the 
measure of a directly accessible physical quantity. Ac- 
cordingly, these different conditions turn the realization 
of an ITP-meter into a highly nontrivial problem. 

Nevertheless it is possible for certain simple cases to 
realize a measurement with an ITP-meter, using for in- 
stance two connected mass transport models as described 
in Sec. IIVBI The systems are similar to the ones shown 
in Fig. [U but now ^2, used as an ITP-meter, is much 
smaller than iSi . We present in Fig. |4] numerical simula- 
tions in which each system Si, {v = 1,2) is homogeneous, 
and obeys the transport rate defined in Eq. ([7T|) with 
fvim) = m^"~^, where r]i ^ r\2- By measuring the den- 
sity pi, of each system, one can determine the value of 
Ai/ using the equations of state, namely A^ = rj^/ Pv In 
practice, one would of course only measure the density 
of the ITP-meter, but here we determine both Ai and A2 
to check the validity of the approach. 

Numerical simulations are run using systems of size 
iVi = 64510 and A^2 = 1024, with equal initial densi- 
ties, pi = P2 = P, and with 771 < 772- Hence the initial 
values of the ITPs are different, namely Ai = 771/p and 
-^2 = Vi/p- At time zero the contact between the two 
systems is switched on, and mass fiows from Si to the 
ITP-meter S2, since Ai < A2. As predicted theoretically 
and confirmed by the numerical simulations, the ITPs of 
the two systems equalize once the steady-state is reached. 
In addition, the value of Ai does not change significantly 
along this process, which is the basic requirement for a 
non-perturbative measurement. Accordingly, the ITP- 
meter indeed measures the value of Ai. Quite impor- 
tantly, this measure is done without knowing the value 
of the parameter rji defining the dynamics of iSi. This 
value was only used to determine Ai from pi in order to 
check the measurement. 

The above example illustrates on a simple model that 
it is in principle possible to measure an ITP using an 
ITP-meter in nonequilibrium systems. Still, in more re- 
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FIG. 4; Measurement process on a system <Si using an ITP- 
meter S2- The values of the ITPs Ai (dashed line) and A2 
(solid line) are plotted versus time. The systems are charac- 
terized by parameters rji — 3, ri2 = 5 (see text). The contact 
is switched on at time t = 0. Initially, pi = p2 = 10, so that 
Ai = 0.3 and A2 — 0.5. Note that Ai almost stays constant 
whereas after a sufficient time A2 converges to Ai. 

alistic situations, finding a suitable definition for the dy- 
namics at the contact that allows for the equalization 
of ITPs turns out to be a major challenge. As seen in 
the preceding section, Eq. (|76p gives a condition for the 
equalization of ITPs between the two systems, and it pro- 
vides useful information to design the contact. Yet, it is 
important to notice that Eq. (ffS)) is a detailed balance 
relation with respect to the stationary distributions of iSi 
and 52- Hence, to satisfy this relation, one requires some 
important information on the gauged system, namely its 
steady state weight fi{q). Such an information is usu- 
ally unavailable in nonequilibrium systems, contrary to 
what happens in equilibrium, where the weights are ei- 
ther uniform (microcanonical ensemble) or given by the 
Boltzmann-Gibbs factor (canonical ensemble). 

2. Measure within subsystems 

An alternative route, that has been exploited recently 
in the context of granular matter [Tol . [26j . consists in 
trying to determine the ITP through an interpretation of 
direct measurements on the system, instead of using an 
auxiliary system (the ITP-meter) [2^ . More precisely, let 
us consider a system with a globally conserved quantity 
Q. Then, from Eq. (jlSp . the variance of the quantity Qat, 
measured over a "mesoscopic" subsystem of size TV, obeys 
the following relation, derived from the grand-canonical 
ensemble: 

{Ql) ^ {Qn? = . (77) 

Let us assume, consistently with the additivity condition 
(HI) , (O , that the variance of Qn is linear in the subsystem 
size N, at least for 1 ^ A'' <|; Ntot, where A^tot is the size 
of the global system. Hence the variance may be written 
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as 



{Qn? = Ng{p) 



with p = {Qn)/N. Then Eq. ^ leads to 

dp 



(78) 



(79) 



Given a reference point (po, Aq), one can then determine 
the equation of state of the system, simply by integrating 
Eq. (l79l) numerically: 



A 



dp 
5(P) 



An 



(80) 



As a first example, we now apply this procedure to a 
mass transport model on a ring with dynamics defined 
by /(to) = (1 + to)2 and v{p) = 1 (see Eq. ((^S])), for 
which the equation of state A(p) cannot be determined 
easily by scaling arguments as in Sec. Ill D[ By measuring 
in a numerical simulation the variance of the mass 
over subsystems of different size iV, we find as expected 
a linear behavior in N already for small sizes (see inset 
of Fig. Efa)), since the grand-canonical distribution is 
fully factorized. The result of the numerical integration 
for the equation of state is in very good agreement with 
the theoretical curve obtained using the grand-canonical 
ensemble (see Fig. [Sfa)), for which p can be determined 
as a function of A, as done for instance in Eq. (jlT)) . 

The fact that this procedure works as well for cases in 
which the steady state does not fully factorize, can be 
seen in another numerical experiment exhibiting a pair- 
factorized steady state. Let the dynamics be defined as 
in Eq. (|34p . with v{p) — 1 and g{in,n) — m + n, that 
is, we assume a — 1, f] = and 7 = 1 in Eq. ((36| . 
Measuring again the variance of the mass Mn over sub- 
systems of size A'', shows that the measurement of the 
ITP A should be performed for higher values of N than 
in the fully factorized case to have access to the linear 
regime (see inset of Fig. EJb)). Note that the results of 
the measurement (see Fig. [SKb)) agree equally well with 
the theoretical curve, which can be obtained in this case 
by the simple scaling argument presented in Sec. HID 21 

In these examples, we obtained for simplicity a ref- 
erence point {po,Xo) using the theoretical equation of 
state. In a more realistic situation where the equation 
of state is unknown, one can estimate as well a reference 
point only through the information obtained from the 
measured variances. Measuring numerically the function 
g{p) for large values of p, one can fit its asymptotic (large 
p) expression with a power law, g{p) ~ Ap^; in the two 
models above, one has S = 2. Then, assuming that A van- 
ishes when /9 — > 00, and that S > 1, the reference value 
Aq corresponding to a given large density po is obtained 
as 



Ao 



dp 



dp 



1 



A{6 



l)Po' 



(81) 
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FIG. 5: Equation of state A as a function of p in two mass 
transport models, obtained by numerical simulations using 
the measurement procedure described in the text (x), com- 
pared with the theoretical prediction p(A) derived in the 
grand-canonical ensemble (solid line); (a) fully factorized 
steady state with /(m) = (1 -1- m)^, (b) pair-factorized case 
with g(m,n) = m + n. Insets: dependence of the variance 
V{Mn) on the subsystem size A'^ (*) for densities p — 1.2 (a) 
and p — 0.5 (b); solid lines are linear fits. 



Note that alternatively, one may also determine a refer- 
ence point from the distribution of Qn, as proposed in 

Accordingly, the present approach provides a rather 
simple way to measure experimentally or numerically 
ITPs in realistic systems for which no information on 
the microscopic probability distribution (like the weight 
factors /(m) for instance) is available. As mentioned 
above, this method has already been used in the context 
of granular material [19, 26]. However, it was thought to 
rely on Edwards' thermodynamic construction which ex- 
plicitely assumes the equiprobability of states compatible 
with the constraints (conserved quantities), whereas this 
assumption is not necessary, as illustrated in the above 
examples -see also ^23]. Indeed, the validity of the ap- 
proach is much more general, making it a convenient way 
to determine practically ITPs in nonequilibrium systems. 



V. CONCLUSIONS 

In this paper, we highlight that the notion of general- 
ized ITPs for steady-state nonequilibrium systems serves 
as a new tool in the study of nonequilibrium phenom- 
ena. Within the range of validity of the additivity condi- 
tion ([5]), ([5]), there is a vast variety of models for which 
our approach can be useful, like models with a factorized 
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steady-state distribution, or one-dimensional systems de- 
scribed by a matrix product ansatz or the transfer ma- 
trix method (at least with finite matrices). The method 
should also be of interest for models which are well ap- 
proximated by a mean field approach. 

Introduced by generalizing the equilibrium concept, for 
which physical intuition is well developed, ITPs benefit 
from a rather clear physical interpretation. They pro- 
vide a convenient way to describe the coexistence of two 
phases, as may be easily illustrated on the case of the con- 
densation transition. An essential issue, which made the 
success of the ITP concept in equilibrium, is the possibil- 
ity to equalize ITPs in different systems put into contact. 
Contrary to equilibrium situations, such an equalization 
is not automatically fulfilled in the nonequilibrium case, 
and the dynamics at the contact turns out to play a ma- 
jor role. We have derived, under simplifying assumptions, 
a coarse-grained detailed balance relation that needs to 
be satisfied by the contact dynamics in order that ITPs 
equalize. Deriving a more general criterion by relaxing 
some of the assumptions made would be interesting to 
gain further insights on this important issue. 

A related difficulty with the approach is that the addi- 
tivity property (jH) , ^ may be hard to test directly. Still, 
measuring the variance of the globally conserved quan- 
tity over subsystems of a homogeneous system may lead 
to an indirect test of this property: indeed, one expects 
that the variance introduced in Eq. (j77p . becomes linear 
in for 1 if the additivity property is satisfied. 

^From the point of view of measurements, it turns out 
that realizing an ITP-meter remains a major challenge, 
mostly due to the difficulties with the choice of the dy- 
namics at the contact. Indeed, one would need to find a 
dynamics for the contact that satisfies the condition for 
equalization of ITPs (that is, the coarse-grained detailed 
balance relation (|76p or a generalization of it), without 
having information on the probability weights of the sys- 
tem on which the measure should be done. Alternatively, 
measurements of ITPs using the variance of the globally 
conserved quantity in subsystems of "mesoscopic" size 
seem to be a promising route. 
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APPENDIX A: INDEPENDENCE OF THE 
CHOICE OF SUBSYSTEMS 



In this appendix, we show that, according to the claim 
made in Sec. Ill Al the value of the ITP Xk does not de- 
pend on the partition chosen. To this aim, let us relate 
the ITP Afc obtained from the subsystems, to a global 



ITP defined on the whole system. Knowing that 

e 

l[dQka'i'{{Qka}\{Qk})^l (Al) 



fc=l 



and neglecting the correction term e^r appearing in 
Eq. ([5]) in the large N limit, one can derive an expres- 
sion for Z({Qk}) depending on the partition functions of 
the subsystems 

r ^ 

Z{{Qk]) = / n ^Q''- Za{{Qka})Zb{{Qk - Qka}) ■ 
k=l 

(A2) 

Consistently with the additivity property (j?]), ([S]), we as- 
sume a general scaling form for the partition functions 
Zn for a large number of sites Ni, 



Zu{{Qku}) = exp<^ A^C.. 



■' aT 



(A3) 



This amounts to assuming that In Z^{{Qku}) is extensive. 
This choice allows us to perform a saddlepoint approxi- 
mation for the integrals in Eq. (|A2[) and thus wc obtain 



Z{{Qk}) = N'Za{{QU)Z,{{Qk-QL})x 



(A4) 



fe=i 



where qk = Qk/N, qu^ = Qkv/N and ql^ = Ql^/N, 
V G {a, 6}. The functions j^, read with n^, — N^/N 



Ia{{qka}) = 



fb{{qka}) 



rib x - 

' 2 ^ 





' qka \ 


^ dqiadq-ja \ 


. "a J 










^ dq^adqja \ 




-q*a){q]a -q*ja) 



(A5) 



Taking the logarithm of the above expression, Eq. (jA4|l . 
leads to 

lnZ({Qfe}) = lnZ,({gL})+lnZ,({Q,.-QL}) 

+ Ine('Zfc) (A6) 

where £,{qk) is given by the finite integrals of Eq. (jA4p . 
We thus obtain for the derivative with respect to Qk, 
using Eq. ®, 



dQk 



dQk 



dQk 



1 d\n^ 



N dqk 
(AT) 

Since £,{qk) is expected to be bounded when N oo 
while the ratio qk — Qk/N is kept constant, the last 
term vanishes in the thermodynamic limit, and we find 



dlnZ{{Qk}) 
dQk 



(A8) 
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where the value of \k stays finite in the Hmit defined 
above due to the scaling form of Z^, in Eq. (|A3p . The 
fact that we can compute the value of Afe directly from 
the whole system shows the independence of the parti- 
tion chosen. This means that the ITPs Xk {k = 1, . . . ,£) 
are global values that can be used to characterize the 
macroscopic state of the whole system. 

APPENDIX B: ADDITIVITY CONDITION FOR 
MATRIX PRODUCT ANSATZ 

The aim of this appendix is to show that, under rea- 
sonable assumptions, a system described by a matrix 
product ansatz with finite matrices fulfills the additiv- 
ity property (O. For the sake of simplicity, we 
only deal with the case of a single conserved quantity, 
but generalizations to several conserved quantities are 
rather straightforward. We also use periodic bound- 
ary conditions, though calculations with open boundaries 
would essentially be similar. Considering a generic lattice 
model as defined in Eq. ((TS]), the conditional distribution 
^(QalQ) reads: 

^{Qa\Q) = ^ TT[Ra{Qa)Rb{Q ~ Qa)] (Bl) 

To proceed further, one needs to introduce the Laplace 
transforms R{s) and M{s) of R{Q) and M{q) respec- 
tively, which leads to R(s) = M{s)^ . Making the fur- 
ther hypothesis that M{s) is invertible, one can find 
a matrix B{s) such that M{s) = exp[_B(s)], yielding 
R{s) — exp[NB{s)]. The matrix B{s) can be decom- 
posed into a sum of two complex- valued matrices |28| : 

B{s) ^ D{s) + L{s) (B2) 

where D{s) is a diagonalizable matrix, and L{s) is a 
nilpotent one, meaning that there exists a positive inte- 
ger p{s) such that L(s)^*^''-' = 0. In addition, the matrices 
D{s) and L{s) commute. Using Eq. (jB2[) . one has 

R{s) = e^^(^) e^^(^) = e^^(^) 7',(iVL(s)) (B3) 

where Vs is a polynomial of degree p{s) — 1, since higher 
order terms in the expansion of exp[A^L(s)] vanish. Ac- 
cordingly, in the large N limit, the dominant contribution 
to R{s) is proportional to Arp('<)-i exp[ND{s)]. Then, as 
D{s) is diagonalizable, there exists a matrix V{s) such 
that D{s) = V^'^{s)A{s)V{s), where A(s) is the diagonal 



matrix. It results that 

e^^(^) =y-i(s)e^^(^)V^(s) (B4) 

so that the dominant contribution to exp[iV_D(s)] is pro- 
portional to exp[N£i{s)], £i{s) being the eigenvalue of 
D{s) with the largest real part. Altogether, R{s) takes 
the following form, to leading order in N: 

R{s) ~ iVP(")-i exp{Nei{s))K{s) (B5) 
where K{s) is a matrix that does not depend on N. Then 
R{Q) is obtained through an inverse Laplace transform: 

j-SQ+ioo 1 

R{Q) - / TT^ NP^'^'' exp[N{£i{s) + sp)]K{s) 

Jsa — ioa ^""^ 

(B6) 

with p — Q/N, and where sq is an arbitrary real num- 
ber, greater than the real part of all the singularities 
of the integrand. Assuming that the equation G(s) = 
d£i/ds -1-/5 = has a solution s*, and that sq can 
be chosen 0(S S , Si saddle-point calculation shows that 
the dominant contribution to R{Q) is proportional to 
exp[N{£i{s*) + s*p)]. The remaining Gaussian integral 
around the saddle (on the imaginary axis) converges since 
G"{s) = £i{s) > 0. This last property can be shown in 
the following way. The logarithm of the grand-canonical 
partition function satisfies InZ(s) « N£i(s), with the 
identification A = s. Then, from Eq. (fT5|) . the second 
derivative of In Z is the variance of Q, which is positive; 
hence f/(s) > 0. Since Z{Q) = Tri?(Q), Z{Q) is also 
proportional to exp[N{£i{s*) + s*p)]. As a result, one 
finds 

lnTT[RaiQa)Rb{Qb)] ~ Nai£U<)+s:Pa) (B7) 

+Nti£\i.sl) + slpt) 
« In Za(Qa)+ In Zfc(Qb) , 

from which the additivity condition ([5]), ([5]) follows. 



APPENDIX C: PAIR FACTORIZED STEADY 
STATE FOR A CONTINUOUS MASS MODEL 

Let us show that the model introduced in Sec. HID 21 
defined through the transport rate given in Eq. 
leads to a pair factorized steady state ((35)) . This can be 
verified with the master equation for this model, which 
reads: 



dPi{m,},t) 
dt 



N 

E 

i=l 



nii+l 



dpip{p\mi^i,mi + p, m,+i - p)P{. 



p,mi 



.0 



/ d^iLp{^\mi_i,mi,mi+i)P{{mi},t) 

._i "'0 
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(CI) 



where P{{mi],t) describes the time evolution of the 
probability of a microstate {rrii}. In the steady 
state the time derivative vanishes. Plugging Eqs. ([34]) 
and ([55)1 into the master equation and dividing by 
Z^^ll^^^g{mi,mi+i) yields: 



AT 

E 

i=l 



5i(mi,TOi+i) g(mj+i,mj+2) 



E 



g{mi_i,m,) g{mi,m,+i) 



(C2) 



Substituting i + 1 by j on the left hand side and knowing 
that g{m,n) — for values of to or n smaller than zero 
one finds the above equality verified. 
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